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Abstract 

We construct the phase diagram of the chiral and diquark condensates at finite temperature and 
density in the 1 + 1 dimensional (2D) two flavor massless Gross Neveu model. The resultant phase 
diagram shows (I) the chiral condensed phase at low temperature and density, (II) the diquark 
condensed phase at low temperature and high density, and (III) the chiral and diquark coexisting 
phase at low temperature and intermediate density. This phase structure is also seen in the 3D 
Gross Neveu model and the 4D Nambu Jona-Lasinio (NJL) model. Thus the phase diagrams of 
the chiral and diquark condensates in the NJL-type models do not change qualitatively in 2D, 3D 
and 4D. 
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I. INTRODUCTION 



Studying the phase structure of Quantum Chromodynamics (QCD) is one of the most 
interesting topics in the field of strong interaction physics. QCD is an asymptotically free 
theory and the interaction between quarks and gluons becomes weak at high temperature 
and density. Then the quarks can be free particles at high temperature and density, and this 
system is called the quark-gluon-plasma. On the other hand, at low temperature and density, 
quarks are confined into hadrons and can not be observed as free particles. Furthermore, 
it is widely believed that the color superconducting phases are realized at low temperature 
and high baryon density. The color superconductivity is the state where the quark-quark 
(diquark) Cooper pairs [2] become correlated, which condense. Since it was indicated that 
color superconducting phases have non-negligible region in the QCD phase diagram^, there 
has been a lot of works on it. For review articles, see Q]. Color superconductivity is 
important because it might be observed in compact stars, such as neutron stars. The density 
of a neutron star is about ten times of that of the normal nuclear. This density corresponds 
to a quark chemical potential \i ~ 500MeV where the deconfinement phase transition may 
happen. Thus it is physically reasonable to expect that color superconducting states are 
realized in neutron stars. 

Investigation of these subjects from QCD is still a challenge to theorists, and we usually 
have to rely on some effective theories of QCD. One of the most important models is the 

n 

Nambu Jona-Lasinio (NJL) model |5| which is a low energy effective theory of QCD. Since 
the NJL model successfully describes the phenomena of the transition from chiral broken 
phase to chiral restored phase, a variety of works has been devoted to the study based on 
it jf|. Recently, it was found that the NJL model also describes the phase transition from 
chiral broken phase to color superconducting phase [zj], and the phase diagram of the NJL 
model shows similar structure with the QCD phase diagram. Thus the NJL model is also 
effective when we consider a color superconductivity. 

The study of the phase structure of the NJL model in lower dimensions (D]_ is also an 
interesting issue, because the model usually becomes simpler in lower dimensions 8j. Indeed, 
the NJL model is renormalizable in 2D and 3D. The NJL- type model in lower dimensions 



is the Gross Neveu (GN) model, which is proposed in 1974 



The GN model shares many 



properties with QCD, such as asymptotic freedom and chiral symmetry breaking in vacuum, 
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and there has been a lot of works on it 
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151 ] . The phase dia gram of quark 



antiquark condensate in the 2D and 3D GN model was obtained in 



and 



131 ] . By using 



the mean-field approximations, the phase structure of the quark-antiquark (qq) and diquark 
(qq) condensates in the 3D GN model with 4-component spinor quarks was derived in [ijj]. 
Although the possibility of spontaneously symmetry breaking is excluded by the Mermin- 
Wagner- Coleman theorem, however, I believe that the study of the qq and qq condensates 
near the Fermi surface is important from the phenomenological point of view. Actually, it 
was shown that the qq and qq condensates could be formed in the 3D GN model by using 
the effective potential analysis [141 ] . The resultant phase diagrams show considerably similar 
structure to the case of the 4D NJL model. This may suggest that the phase structure of the 
quark matter does not change qualitatively in the cases of lower dimensions. To confirm this 
idea, we focus on the phase structure in the 2D GN model here. Recently, within the 2D GN 
model, the qq a nd q q condensates in vacuum (zero temperature and chemical potential) has 



been studied in [15|] , and it has been found that the competition between these condensates 
is sensitive to the ratio of the coupling strengths of the qq and qq condensates. 

In this paper, we study the qq and qq condensates at finite temperature and density in 
the 2D two flavor massless GN model, and obtain the phase diagrams. Then we compare 
the result with those of the 3D GN model and the 4D NJL model. 

The paper is constituted as follows: In Sec. II, we introduce the Lagrangian density of 
the 2D GN model and apply the mean-field approximation. In Sec. Ill, the thermodynamic 
potential is obtained. Then we display the numerical results for the qq and qq condensates 
at zero and finite temperatures in Sec. IV. The phase diagrams of the 2D GN model are 
shown in Sec.V. Finally in Sec. VI, we compare the results of the 2D, 3D GN models and 4D 
NJL model, and discuss the similarities and differences of them. 

II. 2D GROSS NEVEU MODEL 
A. Preliminary 

We consider the following Lagrangian density as the original form of the four fermion 
interaction model, 

£ = qi$q + G s (qq) 2 . (1) 
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Here G,s is the coupling constant of the qq condensate. To incorporate the diquark conden- 
sate in this Lagrangian, we apply the Fiertz transformation to (qq) 2 term and leave the most 
relevant interaction near the Fermi surface [s| . More concretely, the diquark interacting term 
is determined as follows: Let us write the diquark term as (q T Oq) 2 in its most general form, 
where O is an operator in spinor, flavor and color spaces. The Pauli principle, 

q T Oq = O imj = -O ijmi = -q T T q, (2) 

indicates that O is a totally antisymmetric operator. Since the attractive interaction becomes 



the dominant contribution near the Fermi surface js], we select the color antitriplet matrices 
(A2, A5, A7) in color space. The renormalization group arguments 16] force the operator in 
spinor space to be antisymmetric. Then the operator in flavor space has to be antisymmetric. 
Observing the above mentioned consideration, we choose the Lagrangian density of the form 

C = qi$q + G s (qqf + G D ^ (^ 2 A a7 5 g C ) (q C ir 2 \ al 5 q), (3) 

a=2,5,7 

where r 2 is the second Pauli matrix in flavor space and C in q c expresses the charge con- 
jugation. The charge conjugated fields are given by q c = Cq T , q c = q T C, with the charge 
conjugation matrix C . Since a vector in color space can be rotated into one direction 
by a global Si7(3)-color transformation, one can select a color direction to blue without 
loss of generality. This is equivalent to leave \ 2 in Eq.Q, and we have the diquark term 
(qir 2 X 2 -f 5 q c )(q c iT 2 X 2 -f 5 q) . 

B. The model 

Following the reason discussed in the previous subsection, we employ the Lagrangian 
density for the 2D two flavor massless GN model with the diquark condensate term, 

C = qifiq + G s (qq) 2 + G D (qir 2 X 2l 5 q c ){fiT 2 X 2l 5 q). (4) 

In 2D, q = q a k(t,x) is two component spinor fields with the two flavors a = u, d and three 
colors k =r (red), g (green), b (blue). For 7 matrix, we use the following representation, 

n ( 1 \ 1 ( 1 \ 



1 ^ 


l-H 


f 1\ 


U -11 








1 -1 oj 



7 l^o -ij ' 7 l^-i oj = ~ 0, = (5) 

The Lagrangian Eq. (j3J) is invariant under the following transformations: 
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1. Parity V: q(t,x) — > 7°q , (£, — x) 

2. Time reversal T: q(t,x) — > 7°g(— 

3. Charge conjugation C: q c — > — r y 1 q T , q c — > — g T 7 1 

4. Discrete chiral symmetry X^: q — > ^y 5 q. 



The symmetry properties are discussed in the paper 15] in more detail. 



Within the mean-field approximation, the Lagrangian density turns out to be 

1 1 a 2 |A| 2 

£ = qifiq - qaq + - A*(q c iT 2 \ 2l 5 q) + -A(qir 2 \ 2l 5 q c ) - — - i^-. (6) 

Here a and A are the order parameters for the chiral and diquark condensates, 

a = -2G s (qq) and A = 2G D (fiT 2 \ 2l b q) . (7) 

III. THE THERMODYNAMIC POTENTIAL 

A. Derivation of the thermodynamic potential 

The partition function is evaluated by using the standard procedure, 

Z = N' J [dq] [dq] exp j J^dr Jdx(C + fxq lo q) \, (8) 
where /3 — 1/T is the inverse temperature and fi is the quark chemical potential. Introduc- 



tion of the Nambu-Gorkov basis 



id 



* = I 9 I and ¥ = (q f) , (9) 



q c 

enables us to write the partition function in the simple form, 

a 2 |A| 2 



Z —N' exp | — J dr Jdx 



AGs 4G 



D 



X 



fmexJ 1 -^^- 1 )*}. (10) 

"* ^ n,p ' 



The matrix G is defined by 



(r i | (jf-* + W )lfl c *r 2 A 27 5 A n 
*T 2 A 27 5 A* ^-a-^°)l f l c 
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where 1/ and l c are the unit matrices in flavor and color spaces, respectively. 
Then the thermodynamic potential Q = — In Z/f3V becomes 



n(a,\A\) 



<j- 



IAI 



AGs 4G D PV 



In J[<Nf]exp ^^^(/JCT 1 )^ 

L n,p 



where V is the volume of the thermal system. With the help of the formula 



L n )P 



det^^CT 1 ) 



we have 



n(a,\A\) 



a 



IAI 



l-lndet 1 / 2 ^- 1 ) 



4G 5 AG D 13V 
The determinant is calculated by following the same procedure in 

det^CGT 1 ) = [pi - Er] 2 [pl - E A 2 ] 2 [pl - £+ 2 ] [pi - e-*i 



where po = i(2n + l)vrT, (n 



,-2,-1,0, 1,2, •••) and 



E± 2 = (£±//) 2 + |A| 2 , E ± = E±fi,E= J Pl 2 + a 



Then, we obtain 



n(a,\A\) 



a* 



+ 



IAI 



4Gc 4G 



2T EE 



± n 



dpi 

2tt 



ln[/3 2 ( P 2 -E ±2 )]+21n[/3 2 ( P 2 -^ 2 )]|. 



2/2 7-.±2^ 



The frequency summation is performed through using the standard method 



J2HP 2 (pl-E 2 )]=P[E + 2T\n(l + e 



■0E\ 



Thus we finally arrive at 



n(a,\A\) = n (a,\A\) + n T (a,\A\), 



n (a,\A\) 



a 



2 |A| 2 



AG ft AG 



D 



2f d ^[E + El + E A }. 



ln(l + e^ 13 *) + 21n(l + e-" E J; 



Qo is T independent contribution, which is ultraviolet divergent, while T dependent part fiy 
is finite. For the purpose of later use, we rewrite Qq in the integral form in the 2D Euclidean 
momentum space, 



n (a,\A\) = 




The p E terms in the denominators in the arguments of logarithms are inserted so as to drop 
an irrelevant infinite constant. 



B. Renormalized thermodynamic potential 



To eliminate the divergence, we perform the renormalization by introducing the counter 



become 



Z s 
Z D 



Z D 


A| 2 


6 


6 

+- - 


ttc 


IT 


2 


2 

f - 


Tie 


7T 



In a 2 , 
In a 2 , 



(23) 



where e (= 2 — D) is the dimensional regularization parameter and a is an arbitrary renor- 
malization scale. 

Let us introduce the following parameters 



A 



a e 



-7r/(12G s ) 



a e 



-7r/(8Gfl) 



(24) 
(25) 



where o"o (A ) is the gg (gg) condensate at T = = \i in the model in which the gg (gg) 
condensate is absent. By using the above parameters, we finally obtain the renormalized 
thermodynamic potential: 



n r (a, |A|) = n 0r (a, | A|) + Q T (a, | A|) , 
fi 0r (a,|A|) = -(^lna 2 )a 2 -(7rlnA 2 )|A| 2 



4-^ 2 -4iai 2 
pi pi 



Note that the thermodynamic potential Q r has two free parameters (cr , A ). As in 



(26) 



(27) 



14J1, we 



take cr to be the scale of the theory (cr > 0), and choose the ratio A /cr as a parameter. 
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FIG. 1: The chiral( circles) and diquark(triangles) condensates at T = 0. 
IV. THE CHIRAL AND DIQUARK CONDENSATES 

In this section, we show the numerical results for the chiral and diquark condensates at 
zero and finite temperatures. Through minimizing the thermodynamic potential, we study 
the qq and qq condensates. 

In Fig. [Tl we display the results for the qq and qq condensates at T = 0. For A /cr = 1/4 
(panel (a)), we see that the qq condensate is realized for \i = ~ 0.65<t and there dose not 
arise the qq condensate. At \i = 0.7<To, the qq condensate disappears and the qq condensate 
occurs, which indicates the phase transition from qq condensed phase to the qq condensed 
phase. Similar results are obtained for Ao/er = 3/4, 1 where the qq condensate is dominant 
at low fi and the qq condensate dominates at high /x. However in the panel (d) (A /cr = 3/2), 
the qq condensate is dominant for whole /i, and the qq condensate is small compare to the 
qq condensate. Thus the behaviors of the qq and qq condensates are sensitive to the value 
of the ratio Aq/ctq. 

In Fig. El we plot the results at finite temperature for A /ao = 1/4, 1 and \i = 0, 0.4uo. 
For (Aq/oq, fi) = (1/4, 0), the qq condensate starts to decrease at T = 0.25<r and completely 
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FIG. 2: The (^(circles) and (/(/(triangles ) condensates at finite temperature. 

disappears at T = 0.6o"o. For the cases (A /o"o, /i) = (1/4, 0.4<7o) and (1, 0), we observe 
the similar behavior: The qq condensate decreases when T increases and the qq condensate 
does not arise. In the panel (d) (A /cr , /i) = (1, 0.4cr ), the qq condensate exists and it 
decreases with increasing temperature. Thus, the qq and qq condensates decrease when T 
becomes larger, which is the reasonable result because the theory is asymptotically free. We 
have carried out the numerical calculation for the other values of Ao/o"o, and found that the 
behaviors do not alter qualitatively. 

It is to be noted that the qq condensate reduces to zero discontinuously at T = (see 
Fig. [T|). While it decreases continuously at finite temperature for some parameter region. 
These are the signals of the first order and second order phase transitions. The above 
mentioned results indicate that the critical point from first order phase transition to second 
order one appears in the phase diagram, which we will discuss in more detail in the next 
section. 
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FIG. 3: The phase diagrams. 
V. THE PHASE DIAGRAM 

Fig. [3] displays the phase diagrams for the cases A /cr = 1/4, 3/4, 1, 3/2. As seen in the 
diagrams, there appear the following phases: 

C: Chiral condensed phase (a ^ 0, A = 0) 
D: Diquark condensed phase (a = 0, A / 0) 
DB: Double Broken phase (a ^ 0, A ^ 0) 
N: Normal phase (a = 0, A = 0) 

For small A /cr , the C phase is realized at low T and /z, the DB phase occurs at low T and 
intermediate /i and the D phase at low T and high [i. This phase structure bears resemblance 
to the QCD phase diagram. However, when Ao/cxo is large, the C phase no longer appears 
and the DB phase exists at low T and /i, and the D phase does at low T and high /i. Thus, 
as the ratio A /cxo increases, the region of the C phase becomes smaller and those of the DB 
phase and D phase enlarge. 
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The black circles (T c , /x c ) seen in Fig. [3(a) and (b) are the critical points from the first 
order phase transition to the second order one. When T < T c , the transition is of the first 
order and it is of the second order for T > T c . From the figure, we find that the critical 
point moves along the transition line toward n = as the ratio A /cro increases. The critical 
point no longer appears at A /cxq = 1, and the transition from the DB phase to D phase 
is of the first order. The other phase transitions, D — > N and C — > DB, are always of the 
second order. 

VI. CONCLUDING REMARKS 

Through studying the qq and qq condensates at finite temperature and density, we have 
obtained phase diagrams. When the ratio Ao/co is small, there appears the chiral condensed 
(C) phase at low T and /x, the diquark condensed (D) phase at low T and high /i, and the 
double broken (DB) phase at low T and intermediate fi. On the other hand, for large Aq/cxq, 
the DB phase is realized at low T and fi and the D phase at low T and high fi. Thus the 
phase structure drastically changes with respect to the ratio A /cr . Such a behavior of the 
phase diagrams shows close similarity with the ones in the 3D GN model. 

The direct comparison of the present model and the 4D NJL model is difficult since the 
free parameter of these models are different. In the 4D NJL model, the parameter is the ratio 
of two coupling constants Gd/Gs[7], and, in the present model, the parameter is the ratio 
A /<7 . However A is related to Gd through Eq. (l25j) . and it increases as Gp increases. This 
means that A /o"o increases as Go becomes larger. Then, with increasing Gd, the phase 
diagram changes from Fig. [3(a) to (d). After all this, we see that the phase structure bears 
resemblance to the 4D NJL case. 

Having compared the 2D GN model, 3D GN model and 4D NJL model, we see that the 
similar phase structure is realized in all models. This indicates that the phase structure of 
qq and qq condensates in the NJL- type model does not change appreciably in 2D, 3D and 
4D. 
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